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ABSTRACT 


Using radiosonde data from six widely scattered sta- 
tions in the Indian Ocean, an attempt was made to identify 
quasi-periodic wind oscillations during the 1973-1974 
winter monsoon season. The data were analyzed using spec- 
trum analysis and time cross-section analysis techniques. 

In general, the results of the power spectra computed 
show a seven-day peak which seems most significant. A four- 
five day peak and a thirteen-day peak were also found at some 
stations. The cross-spectrum analysis shows little correla- 
tion of the peaks, indicating uncertainty in the significance 
of the spectral results. 

The results of the time cross-section analysis, done 
Over two one-month periods, indicate a possible fluctuation 
time scale of 10-20 days at most stations; and in many cases, 
the fluctuations appear to propagate downward. 

Due to the low density of available stations and the 
poor quality of the data, the results are considered very 
tentative. Further study is suggested, but only after the 


quality of the data can be improved. 
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Te NO ee On 


Synoptic-scale tropical wave disturbances have been 
studied observationally with various techniques. Ordinary 
map analysis has been used in case studies of individual 
Emibems by Riehl (1945, 1967), Palmer (1952), Yanai (1961, 
1963), and others. Such case-study approaches have been | 
highly successful when the data network was adequate. 
Unfortunately, the tropical regions in general suffer from 
a very low density of reporting stations. In such cases, 
a composite approach has been found to be more successful 
than the standard synoptic techniques. Reed and Recker 
(1971) and Williams and Gray (1973) used such an approach 
and found it to be extremely useful in describing the 
average properties of the disturbances. 

A convenient alternative to tne compositing approach 
is the spectrum analysis, which is an entirely objective 
ieemmique. it is important to note that most studies using 
a spectrum analysis approach have been of a heuristic or 
diagnostic nature, rather than purely statistical (Julien, 
1972). Thus, most prior spectrum analysis studies have 
had more in common with the case-study and compositing 
approaches than a purely statistical spectrum analysis. 

The early synoptic studies (Riehl, 1945) showed the 
existence of easterly waves in the Caribbean Ocean. Simi- 


lar waves were also found in the tropical Northwest 








Pacific region by Riehl (1948) and Hubert (1949). These 
waves were found to have a period of ee aan days and 
a wavelength of 1,500-2,000 kilometers. 

Most recent studies in the tropical Pacific Gcean 
have used the spectrum analysis approach, and two groups 
have made significant contributions in this area. One 
group was Yanai and his collaborators at Tokyo University, 
wno used data from mostly the central Pacific stations. 
The second group was Wallace and his collaborators at the 
University of Washington, who used mainly the Western 
faciweac Stations.= Both sroups consistently femmd a wave 
period of four-five days, suggesting that the waves may 
be similar to the easterly waves studied earlier. 

In addition to the four-five day periodicity, Wallace 
and Chang (1969) also found a ten-fifteen day fluctuation 
in the zonal component wavelength on the order of 10,000 
kilometers. 

Spectrum analysis has also been used in other areas. 
Burpee (1972, 1974) did extensive work in the northwestern 
part of Africa during the Northern Hemisphere summer. A 
dominant fluctuation of four-five day periodicity with a 
horizontal wavelength of 2,000-4,000 kilometers was found. 

some observational studies near the Indian Ocean have 
been made but their results generally appear to be quite 
mnconclusive. Rao and Murty (1973), in a spectrum analysis 
study using data from Singapore, Byan Lepas and Saigon, 


found an indication of four-five day fluctuations. However, 
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Peeack Of verticas correlation in this period sandeeaair. 
some doubts to the significance of these PEVicwuia baomor 
Another spectral study ay Saeco and Rao (1974), 
@ver the Indian continent, found periods of eight days 
and three-four days above 200 mb in both components of 
the wind, and a period of ten-twelve days at the levels 
between 900 mb and 500 mb. Malay (1974) did a similar 
study for five Indian Ocean stations and found perio- 
dicities of seven days and fifteen days at several levels. 
Due to the relatively small number of studies under- 
taken and the increased interest in the Indian Ocean, the 
Purpose Of this study is to perform a spectral and time 
cross-section analysis of some newly available radiosonde 
data in an attempt to diagnose the disturbances in this 


region. 
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TIE MOYEN bye 


The study was initially ieee dacivude al ilvawa ie 
Bole tropical Indian Ocean stations. Radiosonde data for 
nine stations were acquired from the National Weather 
Records Center in Ashville. After an examination of the 
Seteimal data, 1t was found only six stations, shown in 
Figure 1, may be considered to have reasonably sufficient 
number of reports to warrant a study. The data covered 
the time period from October 1, 1973 to April 30, 1974, 
which corresponds to the winter monsoon season. These are 
twice-daily reports, which give a total of 424 data points 
for the entire period. 

Due to the many gaps in the data, only five mandatory 
levels from 850 mb to 300 mb were used in this study. In 
addition, only the wind data were found complete enough 
to render the analysis meaningful. 

All the selected six stations were used in two one- 
month time cross-section analysis. However, only four of 
them have sufficiently complete data for the seven-month 


period to warrant spectral analysis. 
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TIT. SPECTRAL ANALYSIS AND RESULTS 


A. POWER SPECTRUM 

The original wind data were placed in zonal (u) and 
meridional (v) components, and the missing data were 
linearly interpolated. Since in several cases as many 
as eight-ten data points were missing in sequence, some 
adverse effect was expected. 

Before calculating the power spectra, the computed u 
and v component time series were filtered to remove the 
effects of fluctuations with periods greater than 25 
days. The filter used was a 98-point high pass Gaussian 
filter designed by Holloway (1958). The standard deviation 
Moers filter is 8.33. The filter reduces the number of 
data points from 424 to 326. The theoretical response 
miumecion of the filter 1s shown-in Figure 2. 

The power spectra and cross spectra were computed at 
the W.R. Church Computer Center of the Naval Postgraduate ‘ 
school. The program used was the BMDO2T of the UCLA 
Biomedical Statistical Program Package. Fifty lags were 
used which correspond to a lag period of 25 days. Spectral 
information for the program was given in intervals of 0.02 
eycles per day (cpd). 

The power spectra were tested for statistical signifi- 
Gance after Mitchell, et al., (1966). Since persistence 


mepresent in nearly all series, it has been found that 
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Theoretical response function for 98-point Gaussian filter with 


m= 8.33 and R(f)=l - exp(-2 n= gé £“), f is frequency. 
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the choice of a Markov "red Nowvse “continuum as cll emeaiauerlan 
continuum is reasonable. Each time series in this analysis 
mas first tested to determine if it fits the null continuum 
or not. This was done by comparing the lag-one serial 
correlation coefficient (r,) With tiaikeor the fircu some 

) or two (5) lags greater than r, to see if the expo- 
nential relation, L5= a and La ry is satisfied. After 
determining that the series does fit the Markov red-noise 
continuum, the null continuum may be constructed by the 
following approximate procedure. Assuming ry is an unbiased 
estimate of the population lag-one correlation coefficient, 
the following eguation may be evaluated for various harmonic 


numbers k, between k=O and k=m: 


Ts po hye 


Lins - 2ry cos #k/m 
In this equation, s is the average of all m+ 1 raw spec- 
tral estimates in the computed spectrum. 

After constructing the null continuum and comparing 
it to the spectrum being studied, deviations of the spec- 
tral estimates were then tested on an "a priori" basis 
to determine if they possess any statistical significance. 
The statistic associated with each spectral estimate is 
the ratio of the magnitude of the spectral estimate to the 
ieecal magnitude of the continuum. Tukey (1950) found this 
to be distributed as oe climes squared ( Xx“) divided be 


degrees of freedom (DOF), where 


DOF = 2N = m/2 . 


mM 
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Here the values of N and m are the number of data points 
and lags, respectively. For this study, the x” / DOF 
value for a one-tailed test of a 95% confidence level 


is 1.68. 


Be RESULTS OF THE POWER SPECTRUM ANALYSIS 

Although the spectral estimates were calculated for 
frequencies up to 1 cpd, only spectral values for frequency 
Hess than 0.26 cpd are shown, because they generally fall 
off rapidly toward higher frequencies. 

The power spectral estimates for the u and v compo- 
nents were plotted and compared with the 95% confidence 
levels as shown in Figures 3 through 6. The confidence 
level is shown as the dashed line in the figures. 

The power spectra for Trinvandrum are shown in Figure 
See For the u component, the 850-mb level (Ug co) shows 
Sremtiicant power in the frequency band of 0.06-0.15 cpd. 
The 700-mb level (Unog) exhibits significant power in the 
Band of 0.05-0.10 cpd. At 500 mb, a single frequency 
0.06 cpd is significant, and at 400 mb a peak occurs at 
Smee epd. Significant peaks are detected in the 0.08- 
aeeee cod band and the 0.17-0.21 cpd band at 300 mb. 

The v-component spectra shows significant power at 
the 850-mb level (Vgcq) in the frequency bands of 0.06- 
meee cpd ands0nd2-0.16%epd, eine Vn 90 spectrum is signifi- 
cant only in the frequency of 0.14 cpd. Neither Veog nor 


Vuog has any significant peaks. In the V300 Spectrum, 


Joes 





-1) 


i 
— 
— wos 
— Wee an 


2 
VARIANCE PER UNIT FREQUENCY INTERVAL (KNOTS“/ 2/50 DAY 





850 ot ee ae eee 850 oer Lj 
mb 0 02.04060810 12 1416.1 820.2224 46 mb O 0204.06081012 1446 18.2022.24.26 
FREQUENCY (CYCLES/DAY) 


Meeure 3. Power spectra for u and v components at Trin- 
vandrum. Dashed line represents 95% confidence level 
for each series, not shown for Viioo* 
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Figure 4. Power spectra for u and v components at Gan. 
Dashed line represents 95% confidence level for each series. 
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Figure 5. Power spectra for u and v components at Singa- 


pore. Dashed line is 95% confidence level for each Serica, 
not shown for Ucog OF Viggo: 
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Figure 6. Power spectra for u and v components at Cocos 
feland. Dashed line represents 95% confidence level. 
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eecienificant band was found in the 0.12-0,16 ced band. 
Over all, the band centered at 0.08 cpd (period 13 days) 
seems to be most significant at this station due to the 


large amount of power in Uso? Yaq9 and some in u 


70 700 


and U3090' 

Subbaramayya and Rao (1974) have calculated wind spectra 
at Trinvandrum. For the u component, they found the spec- 
Geel peaks at eight days (0.125 cpd) and three-four days 
(0.25-0.33 cpd) above 500 mb and ten-thirteen days (0.08- 
0.10 cpd) below 500 mb. For v component, they found spec- 
tral peaks of seven-eight days (0.125-0.14 cpd) and three- 
four days (0.25-0.33 cpd) above 500 mb and seven-eight days 
(0.125-0.14 cpd) and ten-thirteen days (0.08-0.10 cpd) be- 
low 500 mb. Although these periodicity bands seem to fall 
in the vicinity of the present findings, they must be con- 
Sidered inconclusive because of the large variations between 
levels and parameters. 

Figure 4 shows the u and v spectra at Gan. The Ug «0 
and Uno Spectra have no Significant peaks. The next 
spectrum, Ueog? has a lerse peacimethe Owle—O12oucpempame . 
At the 400-mb level, two significant frequency bands ap- 
pear, 0.12-0.16 cpd and 0.20-0.26 cpd. The same bands are 
also significant at the 300-mb level. 

The v component shows some significance at all five 
levels. The significant “peak” shows up in the frequency 
baad of 0.03-0.06 cpd for Vg 50 which is in the filter cut- 


off range and will not be considered. For v, the 


700? 


eps 





Bisniticant frequemcy pand js 0.00-O224+ cpd, andsnor Veo 
moe. O.20-0.25 ,epd.. The Ving Spectrum die) livesiqlahMarorshenre 
mere the bands of 0,.12-0.14 epd and 0.15-0.22 cpd.. At 
V300! it is 0.16-0.24 cpd. Thus, there seem to be two 
common significant frequency bands occurring at Gan--one 
centered near 0.14 cpd (seven days) and the other near 
0.16-0.24 cpd (four-five days). It is interesting to note 
that Malay (1974) has found a significant periodicity of 
seven days at this station for the previous year. 

Singapore's power spectra, shown in Figure 5, have 
almost no significant power in u at all levels. On the 
other hand, the v spectra have significant peaks centered 
ateO.23 epd at 850, 500, and 400-mb levels. The corre- 
sponding common periodicity is about four-five days. 

Figure 6 shows the power spectra for Cocos Island. 
Substantial amount of power in u is seen between 0.08- 
O526 cpd at the three upper levels, 300, 400, and 500 mb. 
The main v peak occurs in the 0.14-0.18 band at 500 and 
850-mb levels, which is also reflected in Ugcg* 

In summary, the significant frequencies seem to vary 
greatly from station to station and sometimes from level 
to level. Thus, no general conclusion can be drawn for 
the entire area; except to say that three period ranges 
centered at the following periods may be significant: 
thirteen days (or 0.08 cpd, for Trinvandrum), seven days 
wee 0.14 cpd for Trinvandrum, Gan and Cocos Island), and 
four-five days (for Gan and Cocos Island). Among them 


the seven day periodicity seems to be most significant. 
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C. CROSS-SPECTRUM ANALYSIS 

Three cross-spectral analyses “eas pe mPormmerl The 
first is the analysis that crosses all levels in the 
vertical at each station for each component. The series 
at 700 mb was always chosen as the base series. The 
purpose of this inter-level cross-spectrum analysis is to 
determine the vertical structure of the fluctuations. 

The second cross-spectrum ahnalysis crosses the u and 
v components at each level at the same station. The 
Wummese Of this analysis is to give an indication of the 
horizontal tilt of the wave, because the phase relation- 
ship between u and v depicts the pattern of streamlines. 

The last type of cross-spectrum analysis performed 
iietars Study was an inter-station analysis. In this 
analysis, the component at one level and one station is 
crossed with the same component at the same ievel but a 
feeterent Station. The results of this type analysis may 
give an indication of the wavelength of the wave and the 
direction of propagation. Obviously, the distance be- 
tween stations becomes an important factor in the inter- 
station analysis. 

| The coherence-squared values and the phase differences 

between two series were given by the cross-spectrum 
analysis. The coherency-squared values obtained were 
fmiest tested for statistical significance utilizing the - 
Peobability points of distribution compiled by Amos and 


Koopmans (1963). Prior to entering the tables, the DOF 


Z4 





were determined from the relationship, 

DOF = 1.25 N/M | 
ie value of N as the number of data points (320) anu 1 
jmseche number of lags (50). The resultant DOF is 6.2 
and the coherency-squared values for the 90% and 95% 
significant levels are 0.280 and 0.348, respectively. 
The phase relationships for the frequency bands with signif- 
icant power spectrum peaks and significant coherency-squared 
values were then used to determine the structure of the 


fine tuations. 


Peek oULTS OF CROSS-SPECTRA 

Generally speaking, only a few cross-spectra were found 
that can be considered as satisfying the significance level. 
For the inter-level cross-spectra, the few significant 
coherency-Squared values all give roughly an in-phase re- 
lationship between different levels. However, the frequencies 
at which these significant correlations occur scatter quite 
widely, and no conclusive results can be drawn from these 
results. The same is true for the inter-parameter and inter- 
station cross-spectra. Their values thus will not be 
presented. It is felt that these poor results may be due 
to the poor quality and especially the large gaps in the 


data. 
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IV. TIME CROSS-SECTION ANALYSIS 


A. TIME CROSS-SECTION ANALYSIS 

Since the spectral analysis provides no conclusive 
results, a time cross-section analysis was performed. The 
time cross-section analysis was done at all six stations 
shown in Figure l. 

Due to large gaps in the wind data, the entire seven- 
month period could not be used. Instead, two one-month 
periods, approximately representative of two monsoom tran- 
Sition seasons, were used for each station. At Trinvandrum, 
Minicoy and Gan, the time cross-section analysis was done 
using October 1973 and April 1974 data; while at Singapore, 
wAaeearta and Coeos Island, the December 1973 °and April 1974 
data were used. The choice of month is based on the fewest 
number of missing data, although it is realized October 
and April are transition months. The u and v components 
were separately analyzed in all cases. 

In this analysis, the average u and v for each cross- 
section were first calculated and are given in Figures 7- 
12. It is seen that most of the monthly averaged wind 
profiles are different from station to station. This may 
be indicative of the complex situation of the monsoon flow 
and may partly explain the different spectral character- 
istics at the different stations. These averaged u and v 


values were then subtracted from the original serics at 
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Figure 7. Average u and v components at Trinvandrum for 


October 1973 and April 1974. 
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Figure 6. Same as Figure 7 except at Minicoy. 
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Figure 9. Same as Figure 7 except at Gan. 
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Figure 10. Average u and v components for December 1973 


and April 1974 at Singapore. 
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Figure 11. Same as Figure 10 except for Djakarta. 
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Figure 12. Same as Figure 10 except for Cocos Island. 
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each station leaving a departure series. Because of the 
uncertainties in the data due to the interpolations and 
partly for the convenience of plotting, the departure 
winds are categorized in nine categories as shown in 


Mapke 1. 


ChTEGORY SPEED (KNOTS) 
0 MEAN 
i _ 07092 2..00 
Z 2.01-4.00 
3 4.01-6.00 
Ly 6 0t—6 700) 
5 6, Ol howeo 
6 10.01-15.00 
7 15.01-20.00 
8 Greater than 20-00 


Table I. Categories of wind perturbations used 
in composite analysis. 


An attempt was also made to investigate any possible 
relationship between the wind fluctuations and convective 
clouds. Cloud amount was determined subjectively from 
NOAA-2 visual images and NIMBUS-5 infrared imagery and 
categorized by the combined brightness of the region be- 
ing investigated, Table II shows the five brightness 
categories. 

CATEGORY BRIGHTNESS 
Clear 
Dark Grey 
Medium Grey 


Light Grey 
Solid White 









Table II. Brightness categories 


215, 





B. RESULTS OF TIME CROSS-SECTION ANALYSIS 

The time cross-sections at Trinvandrum are shown in 
Figures 13-16. The etn u-component analysis (Figure 
13) shows an average periodicity of approximately 13 days, 
which agrees with the spectral findings. MThe fluctuations 
seem to appear oh the upper levels first and spread down- 
ward in time. 

Figure 14 shows the October v-component cross-sections. 
The time scale of the fluctuations appear to be more com- 
placated and the vertical structure is not well defined. 
The appearance of the shorter periodicities may be indi- 
cative of the seven-day spectral peak. 

At first look, the April u-component time cross-section 
(Figure 15) shows no obvious periodicity. However, if 
extrapolation in time is considered, a somewhat longer 
period on the order of 20 days may be assumed. This long 
periodicity may be partly due to the linear interpolation 
of missing data. In this analysis, the fluctuations seem 
to appear at the lower levels first and then spread up- 
ward with time. 

The April v-component cross-section (Figure ANG) Geversulig! 
shows no clear periodicity, only that the time scales seem 
to be shorter than that for u. The vertical structure of 
the fluctuations shows little organization. 

Figures 17-20 show the time cross-sections at Minicoy. 
The October u-component analysis (Figure 17) appears to be 


fairly well organized and shows an apparent periodicity of 
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about 20 days, which again may be due to interpolation. 
The fluctuations appear at the upper levels first and seem 
to spread downward in time. 

The October v-component cross-section (Figure 18) is 
generally disorganized with shorter time scales. No well- 
defined periodicity is obvious. 

The April u-component (Figure 19) and v-component 
(Figure 20) time cross-sections at Minicoy show very little 
organized activity. Neither cross-section shows an obvious 
periodicity. There is no well-defined vertical structure 
clearly visible in either series. 

At Gan, the October u-component analysis (Figure 21) 
shows no clear periodicity. Near the end of the month, 
the fluctuations appear to propagate downward in time. 

In the first few days, high frequencies dominate the fluc- 
miatrons. if one goes back to Figure 13, which shows the 

u fluctuation at Trinvandrum, similar feature can be de- 
tected. Tracking the negative (easterly) u departure for 
day elght at Gan to day nine-ten at Trinvandrum appears 

to suggest a northeastern movement of this particular 
system. However, this feature cannot be detected at 
Minicoy (Figure 17), which is west of Trinvandrum and north 
of Gan, because of the missing data during these days. The 
October v-component analysis (Figure 22) also shows no | 
obvious periodicity. The fluctuations appear to be more 
organized vertically than those found in the u component. 


They also show slight downward propagation in time. 
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The April u-component and v-component cross-sections 
(Figures 23 and 24) show hardly any organization at all, 
except that downward propagation (from 400 mb) is again 
apparent for the negative u anomaly in the last few days. 

Figures 25-28 show the time cross-sections at Singa- 
pore. The fluctuations in the December u-component anal- 
ysis (Figure 25) appear to have an average periodicity 
of approximately 10 days. The vertical structure, in this 
case, appears to propagate upward in time. The December 
v-component (Figure 26) analysis shows about 3.5 waves at 
the lower levels, which indicate a periodicity close to 
seven days. Downward propagation is apparent for this case. 

Tie April u (Figure 27) and v (Figure 28) cross-sections 
at Singapore do not give a clear periodicity. Inter- 
polation of missing data may be responsible for the 
apparent long period at the first half of the period. 

Figures 29-32 show the time cross-sections at Djakarta. 
The December u-component analysis (Figure 29) shows well- 
organized fluctuations with just one wave for the entire 
23-day period, while the April u section (Figure 31) shows 
two waves for the period. The v sections (Figures 30 and 
32) show no well-organized fluctuations. | 

At Cocos Island, all time cross-section analyses 
(Figures 33-36) show well-organized fluctuations. The 
u-component analyses for both months show an average per- 
Mmedicity of about 12-16 days, and so does the v section 


in December. On the other hand, the v-component analysis 
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shows an approximate period of about 24 days. Except for 
the April v component, the fluctuations are nearly in 
phase in the vertical. 

In summary, the results oe the time cross-section 
analyses vary widely from station to station. However, 
in several cases, the u component did show an average 
periodicity of about 10-20 days, where the longer limit 
may be due to interpolation of missing data. On the other 
hand, the v component generally shows shorter scale fluc- 
tuations. In most cases, the fluctuations seem to appear 
at about 400 mb first and then spread downward in time. 

No definite relationship between cloud brightness 
and wind fluctuations can be found. ‘Thus, it seems that 


the two fields are not related. 
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V. CONCLUSIONS 


The power spectra indicate a common peak of about 
seven days at three of the four stations analyzed. fMThis 
periodicity has been observed by Malay (1974) at Guam 
and Subbaramayya and Rao (1974) at Trinvandrum. Other 
possible periodicities found were four-five days at two 
stations (Gan and Cocos Island) and thirteen days at one 
station (Trinvandrum). A period of fifteen-sixteen days 
has been observed in precipitation data over the west 
coast of India (Murakami 1972), which may be related to 
the thirteen-day peak at Trinvandrum found in this study. 
The variation of spectral characteristics between stations 
and levels plus the lack of significant coherence values 
cast considerable doubts to the significance of the spec- 
tral peaks. 

Time cross-section analysis for two one-month periods 
elise provides little solid information. It was found, 
however, that the u component seems to have a longer time 
sealle of fluctuation than the v component. This as con= 
Sistent with observational findings elsewhere, which may 
be a general characteristic of tropical oceanic distur- 
bances. In addition, the fluctuatvons 3m Mesyscsces pp peas 
to begin at upper levels and then propagate downward at 
Bw aay ci mieaiter 

Due to the general poor quality of the data and the 


low density of available reporting stations, very limited 
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results have been found from this study. Further study 
is suggested, but only after the quality of the data can 


be improved. 
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